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I. INTRODUCTION 



In recent years quantum field theoretical progress has been 
achieved towards a reliable description of the dynamics far 
from equilibrium. At the same time, experimental efforts to 
measure the dynamical evolution of complex quantum sys- 
tems in a time-resolved way have multiplied. In particular 
ultracold atomic quantum gases provide a new and unique 
way to induce and observe many-body time evolution far 
from thermal equilibrium. Of primary interest and a funda- 
mental character thereby are experiments setups which al- 
low to induce strong correlations in the system. These in- 
clude, e.g., the time evolution of Bose- and Fermi gases in 
low-dimensional optical lattice traps [2[|, studies of equilibra- 
tion in one-dimensional trapping geometries and the role of 
integrability of the system under consideration JH-Ht]. Re- 
cently advanced experimental system preparation and imag- 
ing techniques promise substantial progress in the study of 
time -resolved many -body dynamics 14j, |a-|9[] . Many interest- 
ing problems in non-equilibrium dynamics have recently been 
studied in other areas like solid-state physics, heavy-ion col- 
lisions, or cosmology. These include, e.g., strongly correlated 
driven situations as transport through quantum dots fioll . the 
dynamics of (quantum) phase transitions Hill , instabilities in 
a quark-gluon plasma IU2I1 . or turbulence in reheating after in- 
flation 111 . 

In this paper we derive dynamic equations for quantum 
fields far from equilibrium, extending our work presented in 
Ref. 01 . We obtain a new time evolution equation for the non- 
equilibrium effective action by use of functional renormaliza- 
tion group (RG) techniques, cf. Refs. Q"j - [l7ll . as well as OH 
for non-equilibrium applications. This exact and closed evolu- 
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FIG. 1 : (Color online) The time evolution of an interacting quantum 
field 0(x, t) (solid line), sketched at some point x, can be determined 
in different ways: By direct integration of an integro-differential 
equation of motion (black arrow), or by the evolution of flow equa- 
tions which, e.g., turn the evolution of a non-interacting field (blue 
dashed line) into that of the interacting one (blue arrows). In this 
article we consider explicitly the causal flow of solutions of the in- 
teracting field equations up to some time r under a change of r (red 
arrow). The flow-equation approach described here can be used to 
determine </>(x, t) in any of the above ways. 



tion equation allows for non-perturbative approximations that 
lead substantially beyond mean-field and quantum Boltzmann 
type approaches: it can be rewritten, in a closed form, as a hi- 
erarchy of dynamical equations for Green functions, and this 
hierarchy admits truncations that neither explicitly nor implic- 
itly rely on small bare couplings or close-to-equilibrium evo- 
lutions. The hierarchy of equations is technically close to 2PI 
and Dyson-Schwinger type evolutions, see I17tl . as well as 
to evolution equations for the effective action as derived in 
lfl9ll . The RG-inspired methods presented here need to be dis- 
tinguished from the time-dependent DMRG techniques build- 
ing on quantum information theory which have been intensely 
studied in the recent past ll20ll . 

As compared to descriptions in frequency space, which are 
commonly employed to obtain, e.g., non-equilibrium informa- 
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tion from spectral densities, real-time approaches profit from 
causality. Causality simplifies the practical solution of com- 
plex integro-differential equations of motion as integrals in 
these equations never involve integrands evaluated in the fu- 
ture, where they are unknown. In the derivation of a dynami- 
cal many-body theory on the basis of functional renormaliza- 
tion group techniques presented in this paper we profit from 
causality in a two-fold way: We first derive the structure of the 
equations which depends on the degree to which interaction- 
induced correlations are taken into account, by means of a 
flow-equation approach. By choosing a cut-off function which 
suppresses fluctuations of the quantum fields at times later 
than the flow parameter r one makes use of causality to derive 
directly non-perturbative integro-differential equations which 
effectively resum higher-order correlations as relevant up to 
the maximum time set by r. Causality then plays a sec- 
ond decisive role in giving a Volterra-like structure to the re- 
sulting equations, allowing them to be solved in an iterative 
way without self-consistency complications between the inte- 
grands and the solution of the differential propagation. 

Our use of the flow equations to derive dynamic equations 
valid up to a given cut-off time represents one possibility. As 
sketched in Fig. Q]it is alternatively possible to derive the dy- 
namical process through a flow from a solution for a non- 
interacting system to the full interacting evolution. 

Functional quantum-field-theoretical techniques which are 
related to the approach presented here and profit from causal- 
ity have been developed and applied intensely in recent years 
O21I - I29I1 . These are based on the two-particle irreducible 
(2PI) effective action i30l - [32ll which is a functional of the 
time-ordered two-point correlation or Green functions like 
(7~\j/t( Xi t)fy(y, t')). Equations of motion for these correla- 
tion functions are derived from the action by means of Hamil- 
ton's variational principle. This procedure requires an explicit 
expression for the action functional which, in turn, is given by 
the series of all 2PI 'vacuum', i.e., external-leg-less diagrams 
which can be written down for a given microscopic model. 

Strong correlations in general require non-perturbative ap- 
proximations. For the 2PI effective action of an 0(Af)- 
symmetric theory an expansion in the inverse of the number 
of field components Af has been studied which in practice is 
obtained by resumming an infinite number of 2PI diagrams 
belonging to a particular class ll2ll - l23i I29I l33ll . In this article 
we show explicitly that one recovers the dynamic equations 
for such a model, in next-to-leading order of the 1/Af ex- 
pansion, without the need of a resummation procedure, from 
our functional flow-equation approach [1]. Beyond this re- 
sult our approach provides a means to derive non-perturbative 
approximations beyond the large-A/" expansion and to check 
the self-consistency of these truncations without the need to 
compute the next order of approximation. Large-A/" approxi- 
mations play an increasing role in experiments with ultracold 
atoms. Strong interest is presently arising in the context of the 
alkaline-earth elements in which the nuclear and electronic 
spins can be almost perfectly decoupled. For this reason, the 
scattering lengths characterizing the interactions in and be- 
tween the Af magnetic hyperfine levels are essentially iden- 
tical, implying a model Hamiltonian with SU(Af) symmetry 



13411 . One has Af — 21+1, which can be as large as 10, 
providing a well-controlled setting for description by as well 
as testing of 1/Af expansions. The Anderson impurity and 
Kondo-lattice models studied in this context have been de- 
veloped in the context of heavy-fermion metals and are thus 
relevant for a large set of applications. 

Our paper is organized as follows: In Sect. [11] we motivate 
our approach and recall basics of the functional-integral de- 
scription of time-dependent quantum field theory. Sect. [Ill] 
gives the implementation of our approach in terms of func- 
tional integrals over a variable closed time (Schwinger- 
Keldysh) contour and details of the derivation of the dynami- 
cal flow equation for the effective action. In Sect.[IV]we derive 
the flow equations for the proper ?i-point functions, n < 4, 
which form the nth moments of the action functional and pro- 
vide the grounds of the derivation of a coupled set of dynamic 
integro-differential equations for the field expectation value 
and the connected two-point correlation function or propaga- 
tor in Sect. [V] In Sect. [VI] we describe in detail an s-channel 
truncation of the resulting dynamic equations and their equiv- 
alence to the equations derived from the 2PI effective action 
in next-to-leading order of a large-A/" expansion. Sect. IVHI 
contains our conclusions, and appendices provide technical 
details of the derivations. 



II. NON-EQUILIBRIUM QUANTUM FIELD THEORY 

A. Motivation 

We consider a system of interacting identical bosons which 
can be described by the Hamiltonian 



H = / d d 



with the one-body Hamiltonian 



H 1B (x) = -—V 2 s + V trap (x), 



(1) 



(2) 



where the field operators obey the usual bosonic commuta- 
tion relations, Vt iap (x) denotes a possible trapping potential, 
and g the coupling constant which, e.g., in d = 3 dimen- 
sions is given in terms of the ,s-wave scattering length a as 
g = Anh 2 a/m. Our approach presented in the following can 
be applied straightforwardly to the case of fermions. 

The dynamics of this many-body system can be described 
in the Schrodinger picture, where the many -body state at time 
t is given by some density matrix po(t). The Schrodinger or 
Liouville equation idtpo(t) = [H, pD(t)]/h, with the Hamil- 
tonian (Q}, determines the dynamical evolution of this state. 
All information about the quantum theory can be encoded in 
the infinite set of correlation or n-point functions 

($t( Xl ) . . . §( Xn ))t = Tr \p D (t)&(xi) ■ ■ ■ §(x n )], (3) 
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with n > 1. For the Hamiltonian (Q]), the equation of motion 
of any particular n-point function involves other correlation 
functions up to the order of n + 2. Hence, the result is an 
infinite system of coupled differential equations for the corre- 
lation functions, in analogy to the BBGKY hierarchy in statis- 
tical physics. This system can in general not be solved exactly 
and requires the choice of a suitable approximation scheme. 

A commonly used leading-order approximation scheme is 
based on rewriting the BBGKY hierarchy in terms of cumu- 
lant or connected correlation functions and neglecting all cu- 
mulants of order three and higher. The resulting set of equa- 
tions for the connected (one- and) two-point functions con- 
stitutes the mean-field approximation which is generically a 
reliable approximation for short-time evolutions and/or small 
couplings: Larger couplings imply shorter evolution times 
over which mean-field approximations are valid. Beyond the 
mean-field level the BBGKY approach yields, in leading or- 
der, to quantum-Boltzmann type approximations. 

In deriving dynamic equations beyond the mean-field ap- 
proximation care is needed in hindsight of conservation laws. 
Crucial quantities like the energy-momentum tensor or the 
particle density-current need to be conserved within the ap- 
proximation of choice as otherwise secularity problems may 
arise and the applicability range of the resulting equations 
be strongly reduced. Functional-integral techniques as ap- 
plied in the following provide means to ensure such conserva- 
tion laws within non-perturbative approximations far beyond 
mean-field and quantum-Boltzmann approximations. 



B. Functional-integral approach 

1. Generating functional 

Before introducing our functional renormalization-group 
approach to nonequilibrium dynamics we would like to re- 
call, for the purpose of making our article sufficiently self- 
contained, some basics about nonequilibrium quantum field 
theory. For more details, cf., e.g., Refs. ll23ll29ll35ll . All infor- 
mation about a non-equilibrium quantum many-body system 
is contained in the generating functional for non-equilibrium 
correlation functions 1 : 



Z[J;p D ] = Tr 



PD(t )Tc exp 



(4) 



Here, 7c denotes time ordering along a closed time path or 
contour C which is discussed in more detail in Sect. Ill B"2l be- 
low. pD{to) is the normalized density matrix describing the 
many -body system at the initial time t , which, in general, 
does not have the equilibrium form ~ cxp{— f3H}. $ 
is the operator of a 2A^-component scalar quantum field, with 
its space-time arguments x = (t,x) = (xq, x). The index a 
enumerates the two components <I> a , corresponding to the real 



and imaginary parts of a single complex Bose field $, or, to 
their independent combinations $ and §\ and it enumerates 
Af different fields describing, e.g., different internal hyperfine 
atomic states. 

Any correlation function of a quantum many-body system 
can be derived from this by functional differentiation with re- 
spect to the classical source field J a (x) and, for a closed sys- 
tem, by subsequently setting J = 0. Of particular interest are 
the connected correlation functions or cumulants which are 
analogously obtained from the Schwinger functional 



W[J;p D ]=-ilnZ[J;p D }. 



(5) 



For example, the field expectation value and the connected 
two-point function or propagator follow as 



4> a {x) = ($ a {x)) = Tr[ PD (t )$ a (x)] 

6W[J;p D ] 



,7=0 



5J a (x) 

Gab(x,y) = {Tc$ a (x)<5> b {y)) - (f> a (x)4> b (y) 
. 5 2 W[J- PD ] 



SJ a (x)SJ b (y) 



(6) 



(7) 



2. Schwinger-Keldysh closed time path 

In Eqs. © and (0, 7c denotes time ord ering along a closed 
time path (CTP) C introduced in Refs. jM®. The time path 
C = C(oo) extends from the initial time to to infinite positive 
times and back to t . We use the short-hand notation 



C(t),x 



dx 



* ,e+ 



dx 



t ,C- 



d d x, 



(8) 



' Here and in the following we use natural units where ti = 1. 



where C ± implies that the times over which it is integrated 
are evaluated on the outward (backward) branch of the CTR 
Note, that in Eqs. ©, @, the infinite real-time contour C can 
be replaced by a contour C(r) extending from to to some finite 
time t > to and back to to- C(r) must contain the times of in- 
terest, i.e., xo and yo, but does not need to extend to any later 
time r > xo and r > £ max = max{x , yo}> m Eqs. © and 
(01, respectively. For any time later than the maximum time 
argument i max of the correlation function considered, the ex- 
ternal currents can be set to zero, J a (x) = 0, xo > t max , such 
that the respective '+' and '— ' contributions on the outward 
and backward branches of C, respectively, cancel in Z, leaving 
C(t 

max ) ■ 

Note also that, if all fields are evaluated at times on the '+' 
branch, the ' — ' branch of C(r) ensures the normalization of 
the generating functional Z[0; pn] = 1, i.e., unity of the trace 
of the density matrix. 

The generating functional Z[J; pn] has a functional inte- 
gral representation which can be found by inserting, to the left 
and right of pn, a complete set of eigenstates of the Heisen- 
berg field operators at the initial time, §{xo — to^x)^^) = 
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Poix)^)' with <Po(.&) = f(to,x): 
Z T [J;p D ] = J [d(p£][dipo](<p£\pn{t )\<Po ) 

x / VipcxpU S T [(p] + [ ip a (x)J a (x) }. 

J<P+,C(t) L L Jc(t).x j ' 

(9) 

We distinguish the generating functionals for different CTP 
endpoints by the subscript r. For r — > oo we regain the origi- 
nal definition |@), Z M = Z. The functional integral 

r ^=/n[ n n 

J v +,C(t) J aS to<a . < r t <a:o<T 

(10) 

sums over all field configurations along the CTP C(r). 5 de- 
notes the classical action, defined as 



C(x), 



(ID 



C(t),x 



where C is the Lagrangian density 2 . We will, for the first 
part of the discussion, not specify this action any further. We 
postpone its definition to Sect.|IV]in which we will also return 
to the non-relativistic Bose gas introduced in Eq. (Q]i. 



3. Initial state 

A general initial density matrix can be parametrized as 

(^\pD{to)\Vo) = Nexp{if e [ip}}, (12) 

with a normalization Af and the functional / expanded in pow- 
ers of the fields: 



/eM = a (0) + E^/ c£U 

n=l >>xi-x n 

n 



(13) 



Here, the coefficients a^].. a „(xi, x n ) vanish identically 
for all times different from to (cf., e.g., Ref. JUt]). 

In many practical cases it is sufficient to specify, at time to, 
only the lowest-order correlation functions. If an initial state 
is fully determined by the mean field 4> a {x) and the propaga- 
tor G a b(x, y), then the initial density matrix, Eq. ( fT2b . can be 
written as a Gaussian in the field <p, i.e., all cr™* 1 with n > 3 
vanish identically. 

We emphasize that non-Gaussian initial density matrices 
can be as easily implemented in the framework we are going 
to discuss in the following. 



to 



FIG. 2: (color online) The closed time path C(r) terminating at the 
time given by the parameter r. At later times, field fluctuations 
summed over in the generating functional do not contribute to Green 
functions the maximum time argument of which is r. 



III. IMPLEMENTATION OF DYNAMICAL FLOW 



A. Variable closed time path 

On the basis of the closed-time-path formulation of the gen- 
erating functional for correlation functions as discussed above 
we will now introduce our real-time flow equation approach 
to non-equilibrium dynamics. The first step will be to define, 
for a particular finite closed time path (CTP) C(r), cf. Eq. (O, 
a regulator function having the form of an additional source in 
the generating functional Z which reduces Z to the finite-CTP 
generating functional Z T . By a Legendre transform of Z T , a 
one-particle irreducible effective action T T [0] is found which 
depends, as Z T , parametrically on the end point r of the time 
contour. This is the starting point for deriving a partial differ- 
ential equation which describes the flow of T T with r. As can 
be seen already from our above discussion of the CTP, this 
flow equation is equivalent to a time-evolution equation for 
the effective action, and therefore for all correlation functions 
derived from it, with time arguments smaller than or equal to 
the CTP's end point. 

Hence, the key idea of our approach to dynamics is to first 
consider the generating functional for Green functions where 
all times are smaller than a maximum time r. The time path 
C(t) closed at t = t, see Fig. [2] leads to the generating func- 
tional in Z T (O. At r = to, this results in a trivial Z to where 
all information is stored in the initial density matrix po{to)- 
From this initial condition, Z T can be computed by means of 
the time evolution d T Z T for all times r > t . 

We will derive this evolution by using ideas from functional 
renormalization group (RG) theory Ill4l4l7ll . To that end we 
note that Z T can be defined in terms of the full generating 
functional Z x = Z in (01 by suppressing the propagation for 
times greater than r. This suppression is achieved by intro- 
ducing an additional two-point source i? T a f,(x, y) which acts 
as a regulator cutting off the functional integration beyond r: 



Z T = cxp < — 



C.x 



y & J a {x 



S 6 
Rr,ab(x,y) 



5 My) 



>Z. (14) 



2 We call S, Eq. (TT}, an action although it is, strictly speaking, a difference 
of two actions, each of which involves a time integration along one part of 
the CTP. 



The regulator function R T is chosen such that it suppresses the 
fields, i.e., 5/6 J a , for all times t > r. This requirement does 
not fix R T in a unique way, and a simple choice is provided 
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fluctuations 




suppressed 


y 


R x (t,f) = 




full path integral 






T 



FIG. 3: (Color online) The cut-off function R T , a b(x, y) in the time 
plane {xo, yo} = {t, t'}, t,t' > to- The function vanishes identi- 
cally except where t = t' > t, x = y, a = b are fulfilled simulta- 
neously where it tends to infinity and therefore implies a suppression 
of fluctuations in the generating functional at times greater than r. 



by 



iRr,ab(x,y) 



oo for xq ~ yo > t,x = y, a = b 

otherwise 

(15) 

see Fig. [3] We emphasize that the above choice fulfills the 
necessary requirements: 

• lim-r^oo R T , ab (x,y) = for all (x ,y ): This implies 
that the full quantum theory is recovered in the r — > oo 
limit, lim T _ ) . 00 Z T [J;p D ] = Z[J\p D ). 

• — iR T is semi-positive definite, ensuring that Z T [J; po] 
does not diverge when derived from a normalized and 
finite Z[J;pd]. Note that the following derivation of 
the flow equations can be achieved without functional 
integrals IU7I1 and that this condition is due to the ap- 
proach presented here. 



—iR r ^ a b(x, y) diverges at r - 
than to, and x = y, a = b. 



> t for x = yo greater 
This guarantees that for 



t = to all fluctuations are suppressed and that Z to re- 
duces to a trace over the initial density matrix po(to), 
see Eq. ©. In this case all ?7-point functions derived 
from Z ta are defined at to only and reduce to the free 
classical ones, or any other physical boundary condi- 
tion we can specify as discussed further below. 

We emphasize that the cut-off R T in (TBl i suppresses any 
time evolution at times greater than r. Correlation functions 
derived from Z T and therefore connected correlation functions 
obtained from 



W T [J;p D ] = -iki Z T [J; p D ] 



(16) 



vanish as soon as at least one of their time arguments is larger 
than t. For example, the connected two-point Green function 



obeys 

G T .ab(x,y) 




for maxjio, yo} > T 
otherwise 



(17) 



Here and in the following, the index r implies that a correla- 
tion function is obtained from a regularized generating func- 
tional, corresponding to C(r). In particular, G T is derived 
from W T in the same way as G = G^ from W, cf. Eq. (|7]i- 

We point out that, in contrast to usual renormalization- 
group flows in the momentum domain, causality prevents the 
influence of times greater than r on the dynamics up to time 
r. Hence, the sharp cut-off < fT3T > is the most physical one in 
the sense that it corresponds to integrating out all fluctuations 
being relevant for the evolution up to a particular time. 

Before we proceed to the derivation of flow equations we 
briefly discuss the possibility of alternative and additional reg- 
ulators. More general regulators R T ,ab{x,y) in the time do- 
main may be advantageous in other cases , e.g., when a cut-off 
is set in the relative-time, i.e., (xo — yo) -direction, see Fig. [3] 
Note that such a case would correspond to choosing a regula- 
tor in the frequency domain lfl8ll which is conceptually closer 
to the traditional RG approach with cut-offs in the momen- 
tum domain. The form of such an alternative temporal cut-off 
could be optimized along the lines described in Refs. lfl6l[T7ll . 

We also remark that additional regulators in the spatial or 
momentum domains may be required to remove divergences 
through renormalization. This can be done along the lines 
discussed in detail in Refs. Ill [Bj Wk ■ 



B. Effective action on finite closed time path 

Summarizing our above discussion, the time evolution of 
the correlation functions of interest can be derived from a flow 
equation for Z T as a function of r. It will be more convenient, 
however, to work with the r-dependent effective action de- 
fined by the Legendre transform 



r T [(t>;p D ] = W T [J;p D }- [ (f> a J a -- [ <p a Rr, a 

JC,a Z JC,ab 



b<Pb- 



(18) 



Here, space-time arguments are suppressed, and we have ex- 
tended the notation © to include multiple integrals, 



/ M ab ..N ab ... -- 

JC{r),ab... JC(r),x a JC(r),x h 

x M ab ( )N ab ...(x a ,x b , ...) , 



(19) 



where M and N are arbitrary functions whose indices a,b, ... 
and space-time arguments x a ,x b , ... are contracted. Possible 
further indices and arguments are not denoted explicitly. 



C. Correlation functions 

In Eq. (fT8l . <p a is the classical field expectation value, 

<f> a {x) = W<#(x)[0;p D ] = SW T [J;pD]/6Ja(x)\j=o, (20) 
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and this relation is to be inverted in order to determine J = 
J[4>; po] used in Eq. ( fT8b . Besides the field <fi we will in par- 
ticular be interested in the r dependent connected Green func- 
tion 



G T ,ab(x,y) = -iW^ b (x,y)[0;p D ] 
. 5 2 W T [J;p D ] 



SJ a (x)SJb{y) 



(21) 



j=o 



The regulator-dependent effective action functional < TT~8T > is 
the generating functional of the regulator-dependent proper or 
one-particle irreducible (1PI) n-point functions 



6 n T T [<l>;p D ] 



i> ai {x\) ■ ■ ■ 5(j)a n (x n ) 



(22) 



(2) 

We finally show that T T ' + R T is the inverse of the propa- 



gator —iG T : From Eqs. (list . d22l we obtain the relations 



T rl( X ) = -Ja{x) - / Rr,ab{.X,y)^) b {y), 
JCy 

5J a (x) 



T rib( X 'y) + R r,ab(x,y) 



(23) 



(24) 



where the functional dependences on <j> and pn have been sup- 
pressed. Taking a further derivative of Eq. (l2Qb gives 



iG T , a b(x,y) 



Hb(y) 
SJ a (x) 



(25) 



Here, J = is taken only after the derivative of <f>. Multiply- 
ing Eqs. (f24l) and (f25]l yields 



F)y ac (x, z) + R T ,ac(x, z) G r ,cb(z, y) 

C,z ' 

S ab 5c(x - y), 



(26) 



where 8q{x — y) = Sc(xo — £/o)<5( x — y) is the 4-dimensional 
delta distribution, in the time variable on the CTP C. 



The flow equation for T T is obtained from Eq. ( fT~8b as 



^T T [<f>;p D ] = -fj-W T [J;p D ] 



i'adrJa 



Cm 



iG T , a t 



CM 



Jb 



4>a(d T RT,ab)4>b 
T Ja 



Cab 

Rr.bc&c 
C.c 



Gr,ab ^rRr,ab 



(28) 



C,ab 



In the first equation above we used Eq. ([25). Terms involving 
functional derivatives of <j> with respect to J, multiplied by 
d T J, cancel each other according to Eq. d23l ). i.e., 



iG-r^at 



C.b 



.(1) 



Jb 



Rr,bc 



C,c 



3 T J a = 0, (29) 



such that the total derivative of T T with respect to r reduces 
to a partial derivative. Using this and Eqs. (|20| | and (|27| | when 
evaluating the derivative of W T gives the second equation in 
d28b . Finally, using Eq. (I261 l. the flow equation d28l i can be 
written as 



d T r T 



f 


1 


ICMb 


Tr ' + R T 



drRr.ab 



(30) 



Eq. d30b represents the central result of Ref . |ljj . This equa- 
tion has the form of Wetterich's equation 111411 for the ex- 
act functional renormalization group flow of the fc-dependent 
effective average action Tk with the momentum cut-off k. 
Hence, it is technically similar to functional flow equations 
used extensively with regulators in momentum and/or fre- 
quency space to describe strongly correlated systems near 
equilibrium Il4l - [l7ll . Its homogenous part relates to standard 
r-dependent renormalization 11711 . and has been studied, e.g., 
in Refs. l39l l40ll . We point out that, in contrast to the usual 
RG flow equations in momentum-frequency space, Eq. (1301 ) 
describes a real-time flow and therefore the causal dynami- 
cal evolution of the quantum many-body system. T T encodes 
the value of correlation functions with time arguments at and 
before r, for a given initial state characterized by r to . 



IV. DYNAMICAL FLOW EQUATIONS FOR 
CORRELATION FUNCTIONS 



D. Dynamical flow equation for the effective action 

The following final step is to derive explicitly the differen- 
tial equations describing the flow of the functionals Z T and 
T T with the cut-off parameter r. We derive the flow equation 
for the Schwinger functional W T by inserting Eq. (fT4l into the 
definition Eq. ( TToT l and using Eqs. d20b and ( TSTI ). 



d T Wr[J; PD } = \ I (G r ,a 
* JC.ab 



b + 4>a<i>b) d T R T .ab- (27) 



A. Model 

The discussion of our flow-equation approach to non- 
equilibrium dynamics has been, so far, very general, in that 
it was relying on only a few assumptions about the underlying 
(bosonic) quantum field theory and a specific choice for the 
regulator matrix R. Extensions to and more specific formu- 
lations for, e.g., fermions and gauge bosons are straightfor- 
ward and shall be postponed to future publications. In the 
following we develop our approach further for the particu- 
lar model of a bosonic many-body system with local quartic 
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self-interactions. We consider a quantum field theory for a 
complex A^-component field tp a {x) (a = (i a ,a), i a = 1,2, 
a = 1, ...,Af; b = (if,, /3), etc.) with quartic interactions 2 , 

~^kj c ^ x )^( x )'Pb(x)(pb(x), (31) 

This model describes, e.g., an ultracold Bose gas of atoms 
with Af hyperfine sublevels whose interaction strength g does 
not depend on the particular hyperfine scattering channel 
('spins') of a pair of atoms. In d = 3 dimensions the cou- 
pling strength in such a system is g = X/Af = Aira/m, a be- 
ing the s-wave scattering length. Hence, the model possesses 
an 0(Af) symmetry under transformations in spin space. The 
free inverse classical propagator of this model reads 



4=0 



5c{x - y)5 aj3 [-iaf alb d Xo - H 1B {x)5ab] 



(32) 



a 2 is the Pauli 2-matrix. We observe that ( l32l is invari- 
ant under transformations iD^ 1 — > UiD^W , with U = 
exp{— ia 2 a}, a£R, such that S T [ip] is invariant under U(l) 
rotations of the complex field <p a = ((fi a + if 2, a)/ V%- 

Hib(x) = — Vx/2m + V(x) is the one-body Hamiltonian 
of a particle of mass m in a possibly time-dependent external 
(trapping) potential V(x). 

Our formulation carries straightforwardly over to, e.g., the 
relativistic Klein-Gordon theory with Xip 4 interactions, in 
which the inverse propagator is fully diagonal and thus in- 
variant under 0(Af) transformations (a, b = 1, . . . , Af): 



iD~b(x, y) = -8 c {x - y)S ab [d 2 Xo -V 2 X + m 2 



B. Dynamical flow equations 



(33) 



The flow equation ( BOb for the full r-dependent effective ac- 
tion T T can be solved exactly only for special models. This is 
not possible for the models introduced in the previous section 
if the spatial dimensionality is greater than zero. The classical 
action for the above models is a polynomial in the fluctuating 
fields cp a , with a single quartic interaction term. One also ex- 
pands the effective action T T in powers of the classical fields 
(\> a and restricts the discussion to the lowest moments T 



(n) 



00 . 
n=Q nl JC, ai ...a n 



3=0 



Here, cf> a is the expansion point, a possibly r-dependent classi- 
cal field. The natural choice is the solution (f> of the equations 

of motion, Tr [<t>] = 0, which extremizes the effective ac- 
tion T r . For this choice, the expansion starts effectively with 



the n — 2 term as the linear term vanishes, and has maximal 
stability. For our model one solution exists where the field 
identically vanishes throughout space-time, 0=0, and one 
or more where </> is nonzero. In the following we will always 
assume that the field <f> is one of the solutions of the equation 
of motion. 

Our goal is to describe the full time-evolution up to some 
time t > t Q , of the functions = T { t n \ for n < 4. In 
particular, we will be interested in the evolution of the con- 
nected two-point function G = ip^ 2 )]" 1 = ipT^] -1 . G 
contains, e.g., the information on the evolution of the single- 
particle density matrix n(x, y, t) = ($^(x, f)$(y, £)) and the 
anomalous density matrix m(x, y, t) = ($(x, i)$(y, i)) (e.g. 
Refs. JMIH). 



1. Flow equation for 

The proper one-point function satisfies the following flow 
equation, which is derived by functional differentiation of 
Eq. <|30]>, 

drTl'M = \ I r ( ? ad M (G T [d T R T ]G T ) cdl (35) 

* JC,cd 

see also Fig. |4] The term in parentheses stands for the reg- 
ularized line. Here and in the following we use the ab- 
breviations ri™2 1 ... a „ = n?o I ...a„(2Cai, . • . ,X a „), G T .ab = 

G T .ab(x a ,x b ), (MN) a b = f c M ac N cb , etc. We point out 

that appears in the equation of motion d23l which deter- 
mines the dynamics of the field cf>. 



(2) 

2. Flow equation for 
The flow equation for the proper two-point function reads 

drT^M = \\ T^abcM {Gr[drRr]G T ) dc 
z JC,cd 



<3) r 

r,acdl 



\G T4e T% f [<t>] 



Cede/ 

x {G T [d T R T ]G T ) fc + P(a,b) 



(36) 



Here and in the following, P(a, b, . . .) stands for a sum over 
the respective previous terms, with all remaining permutations 
of the arguments a,b, . . .. Note that the interchange of, e.g., 
a and b includes an interchange of the field indices and of the 
corresponding space-time arguments. 

To solve the above flow equations, initial correlation func- 
need to be specified which we discuss in 



tions r 



(») 



(34) Sect. lIVB4l below. 



3. Flow equations for coupl ing functions 

The above flow equations need to be complemented with 
corresponding equations for the proper 3 -point function 
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Gz.ab = — © Irlofc-rf = 



jyr^L[*=o]=- 




+ } +- 




will show in the following, they allow to derive, for a partic- 
ular cut-off time t, a set of dynamic equations describing the 
time evolution of the two- and four-point functions up to time 

t = T. 



4. Initial conditions for flow equations 



For the solution of the coupled Eqs. d37b and d3~8l initial 



{ ^J^+P(a,b,c,d)} + ± 

b c 




functions Y 



(2) 
to,ab 



and T 



(4) 

tn, abed 



v v 



FIG. 4: (color online) Diagrammatic representation of the general 



flow equations for Ff 



and T 



(4) 



0], for a 4 - 



theory. Open circles with a r denote d T R T ,ab, solid lines with (blue) 
filled circles are r- and, in general, ^-dependent two-point functions 
G T ,ab = i[Y^ + R T ]2b ■ AH other filled circles denote proper field- 
dependent n-vertices Y "l bcd , n = 3, 4, 6. P implies a sum corre- 
sponding to all permutations of its arguments. Fhe full flow equa- 
tions for Yfl bc [4>] and Y ( f> abcd [<p] are given in App.fA] Fig.© 



, respectively, need to be speci- 
fied at r = to. We point out that, within the truncation scheme 
chosen above, we can insert any set of proper two- and four- 
point functions defined in their time arguments at and only at 
to, as long as we set all n-vertices for n = 1,3, and n > 4 to 
vanish. For higher truncations more n-vertices can be speci- 
fied at to in an analogous manner. 

Fhe truncation scheme discussed here corresponds to a 
Gaussian initial density matrix /Orj(io) since the four-point 
function has an influence on po(t) only for t > to- Here, we 
choose the respective classical proper functions defined by S 
in Eq. (ED, 



T ?abM as wel1 as the 4 -P° int function r^ 6 J4 Fig.g] 
shows a representation of the equations for the r-dependent 



r 



(2) 



ta,ab 



S 



(2) 



ab ' 



(39) 



t„, abed = S abcd = -(VA0(<U<5cd + S ac S bd + 6 ad 5 bc ) 

x Sc(x a -x b )Sc(x b -x c )6c(x c -x d ), (40) 



proper functions Y ( r l \ for n = 1, 2, 4 where the equation for and all hi her n . vert ices vanish. 

( A \ - . n . ^ 



Y^ abcd [(j)] is shown for <j) = 0. In AppendixlAlwe provide the 



full flow equations for Y ( [<j)] and Y ( d [<f)] in diagram- 
matic form. 

Clearly, the equations for these functions involve proper 
functions of order 5 and higher, and so on. For practical com- 
putations, the resulting system of equations needs to be closed 
by truncation at some desired order, which of course is possi- 
ble only at the expense of the exactness of the system. 

We exemplify this here for the case of vanishing field ex- 
pectation value. For = 0, the action (131! implies that 
Tr 3 ' = 0, and thus the flow of Yf' vanishes, see also Figs. [4] 
and|6] Hence, the equation for involves only a term con- 
taining Y^\ 



drT% = - jYZbcAGr \drRr]G T ) dc . (37) 

We need to supplement Eq. ( 137) 1 with the flow equation for 
Y^\ which is illustrated in Fig. [4] (assuming (red) 3-point 
functions to vanish). This system of equations is still exact. 
To close it we truncate by neglecting, in the equation for T^, 
the term involving Y^\ 

dr^rlbcd = ~g jj^tibefGrJg^rldgl^] (G T [d T R T ]G T ) he 

+ P(a,b,c,d). (38) 

Eqs. (f3Tb and ( [381 represent a closed set of integro-differential 
equations for the proper functions up to fourth order. As we 



,(4) 



V. DYNAMIC EQUATIONS 
A. Causality 

Causality implies that no physical observable measurable 
at time t can depend on system properties arising at times 
larger than t. This important physical fact distinguishes quan- 
tum field theory on the real time axis and the flow-equation 
approach from their formulation for equilibrium situations. 
Causality manifests itself, in the present theory, for the sharp 
cut-off regulator defined in Eq. ( fl~5T > and the temporally local 
models introduced in Sect. IIV Al in the following properties 
of the correlation functions. 



G Tt ab — i 

(G T [d T R T ]G T ) ab = 



r (2) 

x r 

r (2) 



(T-t a )8{T-t b ), (41) 



d T [6(T-t a )9(T-t b )]. 



(42) 



For easier readability we shall use the notation t a = x a fi, etc. 
The ^-functions are defined such that 0{T—t) is zero for t > r, 
otherwise it is 1 . Relations similar to Eq. ( HTI) can be derived 
straightforwardly for all connected n-point Greens functions, 
see AppendixlBlfor the details. 

We have chosen the sharp cut-off regulator ( TOl ) in accor- 
dance with causility. As discussed above, the flow equation 
for the generating functional Z T , and therefore the flow equa- 
tion for the effective action Y T are time evolution equations for 
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the respective functionals. Hence, the dynamics of the proper 

(n) 

vertices If is readily given by the flow equations derived in 
the previous section. 

At this point different avenues can be taken to solve these 
equations, see Fig. Q] For the regulator ( fT31 ), a convenient 
procedure is to integrate analytically the flow equations for 

(n) 

the Tt , from initial time to to some present time t > to- 
The resulting equations have the form of integro-differential 
equations for the n-point functions, representing initial-value 
problems. The form of these equations is familiar from other 
approaches discussed in the literature so far. We adopt this 
approach in the present article. 



B. Dynamic equation for the field expectation value 

We begin by deriving the equation of motion for the field 
expectation value 4> a {x a ) defined by the solution of 



r (i) 



(43) 



Here we have used that the flow vanishes for r > t a = x a> o, 
see JB3t . Note also, that then (l43l coincides with d23l at van- 
ishing external current. In order to derive the dependence of 
the one-point function 1^ a [</>] on its initial condition [</>] 
and the time evolution of the correlation, we decompose the 
effective action Ft [<f>] as 

t 

T t [4>]=T t0 [<f>] + J drdrTM, (44) 

to 

and calculate its functional derivative with respect </> c : 



64>c(x c ) 8(j>c{x c ) 



At 



84> c {x c ) 



Cab 



dr 



J ab 



,(3) 



MGr[d T Rr}G T ) ba . 



la 



Cab 



(45) 



Using the identities (HTI) and d42l) we can carry out the integra- 
tion over the cut-off time r. We obtain the dynamic equation 
for <f) c by inserting Eq. d43T > into d43l ): 



i 



^ b ,caM G Tab ,M = 0- (46) 

z ■IC(t a ),ab 



Here, we have introduced the notation 

r ab = max{t a ,t b } 



(47) 



as a shorthand for the maximum of time arguments of the two- 
point function. We finally insert again the classical action ( l3TT i 



as initial condition to the integrated flow equation (l46b and 
arrive at the integro-differential equation of motion 



C(t a ),d 
1 



iD cd ( l } dj ~ -^j(l>d{Xc)4>d{.X c )(j) c {x c ) 



T { ^ C M G TahM m. (48) 



2 JC(t c ),ab 



All time arguments of correlation functions involved in this 
equation are smaller than the present time i max = t c . This 
can be shown by considering the integral equation for T^ b 
provided in Appendix IB 41 Hence, Eq. (08]) is causal and can 
be used with the set of equations derived in the following sec- 
tions to solve the nonequilibrium time evolution of a system 
iteratively. 

The time evolution of the field expectation value cf) c is also 
encoded in the dynamical flow of the one-point function d35l ). 
We take the r-derivative of the equation of motion (l43l and 
arrive at 



JCd 

Solving d49l ) for d T (f> results in 



(49) 



dAc = i / G T , cd d T r^ d [ct>] 

JC(r),d 

= ~\ I G TtCd T^ dab (G T [d T R T ]G T ) ba . (50) 

Z JC(T),abd 

The r-evolution of cf> is complementary to the time evolution 
of (f> for times smaller than r. Eq. (T5Qb completes the set of 
flow equations for n-vertices that is evaluated on the equation 
of motion. 



C. Dynamic equations for propagator and 4-vertex 

1. Dynamic equation for the propagator 



In this section we use the dynamical flow equation ( I36I I to 
derive a integro-differential dynamic evolution equation for 
the proper 2-point function. Integrating Eq. (f36l > (recall that 
C = C(oo)) and employing Eqs. d4"TT ) and d42b yields 



,(2) 
T,ab 



to 



to 

t 

\ I ^rlbcd (Gr[d T R T ]G T ) dc 

L J JC(oo),cd 



p(3) si r yo> 

L T,acd °>,de 1 T ^ be f 



(3) 



C(oc),cdef 



(G T [d T R T ]G T ) fc + P(a,b) 
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- f r (4) r , 

9 / L T cd , abed '-'■Tad, dc 

2 JC{t),cd 



,(3) 



2 Jc(t), 



cdef 



T adc f,acd ^Tcdef, 



x r (3) g 



(51) 



Here, we have used the notation ( l47l ) and, analogously, 

Tabcd = max{i a , t bl t Cl t d } (52) 

as a shorthand for the maxima of time arguments of the re- 
spective two- and four-point functions. d5"ll is the dynamic 
equation for the two-point function G ab for the evolution from 
to to t. We will explain this in more detail in the following. 

In the case of t being smaller than the maximum d4Tb of t a 
and t b the right hand side of Eq. dSTb must be zero because of 
the causality property jB lb . Due to relation dB31 >. we can set 
the cut-off time to any value t > t a ,t b . 

From the above integrated flow equation we determine the 
integro-differential dynamic equation for the propagator by 
contracting it with [v[ 2 ']^: 



C(oo) ,a' 



.(2) 
tQ,aa' 



,(2) 



5 e {x a - x b )5 ab 



C(oo),a' JC(t),cd 



,(4) 



r c rf,aa' cd T cd 



G 7 



, d c 



,(2) 



C(oo),a' JC(t),cdef 



r (3) r 



x r (3) r 

X 1 T cd ef,a'ef U T cd ef ,/c 



,(2) 



a' b 



(53) 



Our focus lies here on the microscopic theory defined by the 
classical action Eq. ( l3Tb which is local in time. Hence, the 

(2) i 

inverse classical propagator F£ ab = iD ab is diagonal in its 
time arguments. For a theory of this kind, the proper vertices 
satisfy Eq. ( IB II ). and thus the integration over t a i on the right 
hand side of Eq. d53l is restricted to the time path C(t). Choos- 
ing t > t b and using identity (PiTl we can replace [rp^]~,^ by 
—iGta'b m the above equation. 

Inserting, moreover, the initial condition d39l ) we finally ob- 
tain the dynamic equation for the propagator in the form: 



C(t),a' 



iD aa>Gt,a'b = iSc(x a ~ X b )S ab 



f 

O / 1 T cd ,aa'cd y - jT <:d,dc^ r t,a l b 

z JC(t),a'cd 



'cd 
,(3) 



C(t),a'cdef 



T^efMcd^T^e-f.de 



x r 



(3) f~i (~t 



(54) 



The time derivatives in this integro-differential equation are 
embodied in r^ o6 = iG^ ab , see Eqs. © and d33l l. 



2. Integral equation for the vertex 

The dynamic equation (154-b needs to be complemented with 
equations for the three- and four-point vertices. In the fol- 
lowing we will again consider first the special case that the 
field expectation value <p vanishes, see Sect. lIV B] and present 
the result for the integrated equation for T^ abcd for this sim- 
plified case. In order to obtain a closed set of equations we 
furthermore truncate the equation for the vertex by dropping 
the diagram involving the six-point function, see Fig. [4] 

The derivation proceeds in the same way as in Eq. dBll l for 
the proper two-point function and yields 



r 



(4) 



ibed 



^(4) 



C(t),efgh 



e fgh ,abefGT fg ,fg 



X T r^,cd g h G r eh ,he + (a O c) + (a O d). (55) 

The brackets denote terms with the given indices swapped. 

3. Implications from causality 

Let us, at this point, come back to the consequences of the 
causality properties found in Sect. IV Al for the dynamic equa- 
tion (|54l and ([55). Without loss of generality we can assume 
t a > h since in the reverse case, the propagator can be ob- 
tained by symmetry G T , ab {x a , x b ) = G rM (x b , x a ). 

We thus can choose t a to coincide with the present time 
i max at which Eq. ( 1541 determines the further time evolution 
of the propagator, and set t = r ab = t max . Then, all time 
arguments appearing in the integrands on the right hand side 
of Eq. d54"b are evaluated at times < i max - 

The truncated dynamic equation (1541 then reads 



C(r ab ).a' 
1 

2 



iD aa l( ^r ab M'b = iSc(x a - X b )8 ab 



^ { r] d Ma'cd G ra d .dcG Tab ,a'b- (56) 



C(r a b) ,a' cd 



Note that neither in Eq. d55l nor in d56l . the flow parameter is 
any longer a free parameter. It is fixed by the time arguments 
of the correlation functions. In particular, the flow parameter 
of the propagator is everywhere set to the maximum of its 
time arguments and could be dropped, such that G T = G is 
recovered as the full Green function. 

The time arguments of the proper four-point function ap- 
pearing in Eq. ( l56l and determined by Eq. (1551 are also at 
most as big as the maximum time r ab = t max . In the integral 
equation ( l55l for the vertex the cutoff parameter is evaluated 
at the maximum time T e f g h occuring in the loop. Moreover, 
taking into account the relations dB2b and (IB3b following from 
causality, the cutoff parameters r e f g h on the r.h.s. of Eq. 031 
can be replaced by either to or the maximum time of the re- 
spective four-point function. 

This leads to a fully causal structure of the full set of trun- 
cated dynamic equations exhibiting the full strength of our 
approach: The dynamic equations for the evolution from to 
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to tmax can be used to solve the time evolution of the system 
iteratively. At the maximum time t no information is needed 
to propagate the evolution than what has been computed in 
the past. In the following chapter we explicitly show that it 
is this causal structure of our equations which allows to red- 
erive from them, in an s-channel truncation, the nonperturba- 
tive dynamic equations known from the 2PI effective action in 
next-to-leading order of a large Af approximation. 



zero during the entire time evolution. This property of the 
flow equations is a consequence of the Z 2 symmetry of the 
microscopic theory. The effective action remains in the sym- 
metric phase, i.e., = 0, during the flow if it is initially in 
the symmetric phase. Thus and all odd proper ?i-point func- 
tions remain zero. The following discussion will first focus on 
the symmetric phase and later be extended to the symmetry- 
broken phase 7^ 0. 



VI. NONPERTURBATIVE s-CHANNEL TRUNCATION 

We are now in the position to put the formalism to work. 
In this section we evaluate the system of equations for 0, G, 
and within a non-perturbative s-channel approximation. 
This allows us to derive a closed system of dynamic equations 
which are then solved numerically with presently available 
hard- and software technology and which goes substantially 
beyond Hartree-Fock mean-field and Quantum-Boltzmann ki- 
netic equations. We show, in particular, that the form of 
the equations for the propagator corresponds to that of the 
Kadanoff-Baym dynamic equations and that the s-channel ap- 
proximation of the self-energy functional is identical to that 
derived from the 2PI effective action in next-to-leading order 
(NLO) of a 1 /TV expansion JHHH. 

As above we consider the complex scalar model defined 
by the classical action ( TUT i and discussed in more detail in 
Sect. IIV Al Most of the discussion presented here is inde- 
pendent of the particular choice, both of Af and D~, . We 
will focus on the case of Af complex fields describing an ul- 
tracold Bose gas of atoms in Af internal states with identical 
s-wave contact interactions. Nonetheless, the derivations pre- 
sented will apply equally well to the relativistic scalar 0(Af)- 
symmetric theory the free inverse propagator of which is given 
inEq. ( f33l >. Hence, our results concerning the 1/A^-expansion 
of the 2PI effective action are directly relevant for the dynam- 
ical theories presented in the work cited above. 

As before we consider Gaussian initial conditions, i.e., only 
the proper functions up to order n = 4 are assumed to be 
non-zero at the initial time to of the flow, see the discussion at 
the end of Sect. HVBl The truncation scheme we use was de- 
scribed in detail in the previous section and takes into account 
the flow of the proper functions up to n = 4. 

Before proceeding we consider the flow of the three-point 
vertex, see Fig. [6] which is proportional to the three-point ver- 
tex itself provided one neglects, as in our truncation, the five- 
point function. For our model defined by the classical action 
OH the bare three-point function is given by 



S 



(3)[ 



A 



dc(x a - x b )Sc{x b - x c ) 



X [S a b4>c(x c ) + Sac4>b{xb) + 8bc4>a{x a )\ (57) 

and is therefore proportional to 0. Thus, if we start with a field 



expectation value = and with r^ a{lc [0] = S^l' c [(j)}, the 
three-point function is identical to zero for all cut-off times. 
Since the time derivative of is, according to Eq. d48l ), pro- 
portional to the three-point function and 0, the field remains 
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A. Dynamics in the symmetric phase 

1. Dynamic equations 

In the s-channel approximation introduced above we are 
left with the closed set of flow equations for the proper two- 
and four-point functions, Eqs. (fJH and (|55j, respectively. 



,(2) 
T,ab 



,(4) 
r.abcd 



1 / ^Tc d ,abcd,GT cd Ac, (58) 

1 JC(t),cd 



G 

f„ 9 L,^ f -Tef 9 h,abef T fgJ9 r^ Sah ,cdgh 
z JC{t),efgh 

G T<lh ,he + (a^c) + (a^d). (59) 



,(4) 



The resulting dynamic equation for the connected two-point 
function G is given in Eq. 



2. Mean-field approximation 

We start by considering the mean-field or Hartree-Fock 
(HF) approximation. The mean-field limit implies that the 
evolution equations are local differential equations, i.e., the 
propagation kernel of any correlation function is local in the 
time and spatial variables. We obtain the mean-field approx- 
imation by neglecting the flow of the four-point function. 
Hence, keeps, for all r, its initial form set in Eq. d40b . 
For this truncation, the dynamic equation for the propaga- 
tor reduces to the well-known time-dependent HF equations 



i D ac{ x a,x c )G Tab , cb (x c ,x b ) = i5 ab 5 c (x a - x b ) 



C(T ab ),C 

A 



G ta ,dd(x a , X a ) 8 ac + G ta ,ac(x 



(60) 



These are differential equations in the first time argument of 
the two-time Green function G. They can, however, be re- 
duced to differential equations for the time- diag onal functions 
G ab (x,y) with xo = yo, see Refs. 11261 12911 . i.e., the off- 
diagonal elements with .To 7^ yo are irrelevant for the evo- 
lution on the diagonal. Usually, the HF equations are derived 
by the use of the equation of motion for the two-point Green 
function, under the assumption that the four-point function is 
a product of two two-point correlators 14211 . This Gaussian 
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decorrelation assumption is at the heart of the above neglec- 
tion of the flow of T^ 4 \ We emphasize that due to the missing 
two-particle (four-point) correlations, the life time of small 
excitations of the system is infinite, and therefore the HF equa- 
tions do not account for scattering processes which could lead 
to equilibration. 



3. 



-Channel approximation 



It is necessary to consider the flow of the four-point func- 
tion to derive dynamic equations that contain real collision 
effects of the particles. In this section we shall consider an s- 
channel approximation to the truncated flow equation of the 
four-point vertex I^ 4 '. We will show that this approxima- 
tion is equivalent to the widely used NLO l/Af approxima- 
tion of the 2PI effective action l2lH23l l29tl . We emphasize 
that, however, the derivation of these equations within our ap- 
proach does not require a perturbative resummation of an in- 
finite number of diagrams as it does within the 2PI approach. 
Moreover, rederiving the dynamic equations in NLO of the 
2PI l/Af does not only reproduce an already known result but 
offers possibilities to obtain new insight into the nature of the 
approximation. 

Within the present approach the s-channel approximation 
is based on the assumptions 



p(">4) ^ r> 
x T,ai...a„ u ! 



(~\ — p wi" 1 ') <r p r( m ) 

I 1 r s ) -L T,a!...a m r s 1 T ,ai...a m i 



(61) 
(62) 



for to = 2, 4. Here, P s projects the proper two- and four- 
point functions onto the s-channel structure of the two- and 
four-point vertices, respectively. The flow-equation approach 
provides the possibility of a self-consistency check of the as- 
sumptions doTt and d62l l. This check involves calculating, us- 



ing Eq. ([59j, the terms (1 - P s ) F^ b and (1 - P s ) T 



,(4) 



as well as V 



(0) 



" T,abcd' 

t abcdef tnat aie generated by the s-channel 
proper two- and four-point functions. The above assumptions 
are self-consistent if the generated functions are sufficiently 
smaller than the corresponding s-channel terms. 

We now define the s-channel approximation. It takes into 
account, in addition to the Hartree-Fock terms, the flow of one 
scattering channel in the equation for the four-point function. 
Labeling the field indices and space-time arguments of the 
four-point function T^ 4 ) in a unique way, the function has s- 
channel tensor structure if only elements with pairwise equal 
field indices and space-time arguments are nonzero. The pos- 
sible combinations of pairs are usually referred to as the s-, t- 
and w-channels. Any nonzero four-point function with three 
or more different indices or arguments is excluded from the 
s-, t-, and M-channel structures. The four-point vertex in ei- 
ther of these channel projections is no longer fully symmetric 
under the interchange of any two of its indices and arguments. 



Note that the classical four-vertex T 



(4) 

tQ,abcd 



sltld which is 



given for our model in Eq. (|40l ) consists of an s-, a t-, and a u- 
channel term, often called direct and exchange terms. While 
all three channels are taken into account on the mean-field 



5*4 




(to U= %...% ) 



(t d = t ...mxt(t a ,t b ) ) 



FIG. 5: (color online) The upper equation is the s-channel projection 
of Eq. J55 b given in Eq. l !65t . The second equation defines the re- 
summed vertex appearing in the Dyson-Schwinger equation derived 
from the NLO l/Af approximation of the 2PI effective action l2lll . 
The two definitions are identical in every order of a perturbative ex- 
pansion (see text). Dashed lines denote the s-channel part of the 
bare vertex ab , see text. All other symbols correspond to those in 
Fig. [4] Letters on internal lines indicate summation over field indices 
and integration over space and time (along the CTP from to to the 
time given in parentheses and back). 



level, an s-channel truncation is chosen for the flow of the 
four-point function Eq. ([59}, by inserting 

X 5 C (x a - X b )6c{xb - X c )Sc(x c - Xd), (63) 

into Eq. d59l ). The integral equation for the s-channel vertex 

^ T.abcd(^ a ' *^ c ' ^d) — &ab$cd 

x 8 c {x a - x b )5 c {x c - x d )Tf )s {x a ,x c ) (64) 

reduces to 

^ {Xa j Xb) 



t _ i 

to ~ 2 



r|- e ^ (^a;*£e) (*T eg ,eg i$e j ^g) 

C(t),eg 

X TlffiXg^b) G Tcg , ge (x g ,x e ). (65) 



for a vertex function t[ 4 ^ s (x a , x b ) which is a scalar in field- 
index space. This integral equation is depicted in Fig. [51 upper 
line. 

Before we discuss the structure of this equation further we 
provide the s-channel projection of the dynamic equation ( 
for the two-point function G T . a b(x a ,Xb)' 



*f ac {Xa,Xc} ^-*ac( x a) ^T ab ,cb{x c , Xb) 



C(r ab ),c 



i5 ab 8c(x a - X b ) - i-^j 



C{r ab ),c 

X Gr ac ,ac{Xai X c ) G Tcb-C b( < X c: Xb) 

where 



(66) 



G ta ,cc{x a , X a ) Sab + G ta 



A 

2Af 



(67) 
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is the HF mean-field contribution to the self-energy shift, cf. 
Eq. ( 160b . and 



i t (x a ,xb) = -ijr^fij.ij) 



(68) 



Equation d66l l together with the definitions d67b and ( l68l ). 
and with the integral equation do3T l for the coupling func- 
tion describe the time evolution of the propagator G in the 
s-channel truncation of Eqs. (l58l l and d59l ). These equations 
depend only on the initial conditions D^ 1 and and can be 
expressed solely in G a b(x a7 Xb) = G Tab>a b(x a , Xb)- Nonethe- 
less, the coupling function (x a , Xb) depends, besides the 
arguments x a , Xb and the indices a, b, on the cut-off parameter 
t which takes independent values, e.g., in Eq. ( 1651 ). 

In Appendix [C] we show that the coupling function and 
therefore the above dynamic equations can be rewritten such 
that they get rid of this apparent additional degree of freedom 
and that they are equivalent to the dynamic equations derived 
from the 2PI effective action in NLO of a 1/Af expansion, 
cf. Ref. I2H1 and Fig. [5] lower line. The resulting coupling 
function is defined in Eqs. ( lC9l >, dC16t , and ( IC17I ). The set of 
2PI dynamic equations in NLO 1/Af approximation has been 
used extensively to investigate nonequilibrium dynamics and, 
expressed for F and p, can be found, e.g., in Refs. ll2ll - l23Tl29tl . 

We emphasize that the proof builds crucially on the causal 
structure of the dynamic equations derived in our flow- 
equation approach. This manifests the strength of the flow 
equation approach originating in causality as already pointed 
out in Sect. IV C 31 above: No information is needed for the 
propagation of the dynamic equations at time t which has not 
been computed during the previous time evolution. This ren- 
ders the evaluation of the integral equations technically feasi- 
ble. 

Complementing the formal proof in Appendix [C] we show 
in the following that the above identity can be inferred in a 
comparatively easy way from the topology of the different 
terms in the flow equations for the two- and four-point func- 
tions: Consider the untruncated set of equations as displayed 
in Fig. |U First, non-s-channel contributions do not generate 
bubble-chains of the form shown in Fig. [5] Second, the term 
involving the six-point function to the flow of T y T ' does not 
give rise to bubble-chains if inserted recursively into the first 

( A\ 

diagram on the right-hand side of the flow equation for T T . 
To see this one has to generate the six-point function from the 
s-channel approximation for the four-point function. We em- 
phasize that beyond the s-channel approximation the six-point 
term in the flow of the four-point function does create bubble 
diagrams. In turn, by dropping the second diagram with the 
six-point function and using the s-channel truncation, the iter- 
ated flow equation generates only bubble-chain diagrams with 
full propagators as lines. Hence, a r-integration of this set of 
flow equations leads to dynamic equations which include all 
bubble-chain contributions containing only bare vertices and 
full propagators and are therefore equal to those obtained to 
NLO in a 1/Af expansion of the 2PI effective action 1I21I I22I1 . 
We emphasize that the above topological arguments are gen- 
erally valid when comparing resummation schemes inherent 



in truncated RG equations of the type of Eq. (f3Ql >. with those 
obtained from 2PI effective actions. This applies, e.g., to equi- 
librium flows lfl5l - [l7ll and thermal flows ifla l43ll . For a com- 
parison with 2PI results see Ref. 14311 . for the interrelation of 
2PI methods and RG flows Ref. 11711 . Note that decomposing 
the six- and higher n-point function according to their flow 
equation into diagrams containing only bare vertices and full 
propagators yields a loop expansion which, truncated at a par- 
ticular order, is inequivalent to the expansion obtained without 
these higher-order vertices. This is analogous to the general 
inequivalence of loop expansions derived from nPI effective 
actions with different n 14411 . 



B. Dynamics in the symmetry-broken phase, (f> / 

1. Dynamic equations 

The complete flow equations for the one- and two-point 
functions are shown in Fig. [4] In the case of spontaneously 
broken U(l) symmetry, i.e., for a nonvanishing field expec- 
tation value 4>, the three-point vertex Tf 3 ' is in general non- 
zero. Hence, we need to take into account the flow equations 
given in Appendix lAl see Fig. [6] As before, we will consider 
only truncations, where we keep proper n -point function up to 
n = 4, such that all diagrams will be omitted which contain a 
proper five- or six-point function. This results in the dynamic 
equations for ri 3 ^ and given in App. lB 4| see Eqs. (IB 101 ) 
and (IB lib , respectively. The dynamic equations for the propa- 
gator and the field expectation value were defined in Eqs. 
and (1481) . respectively. 



2. Mean-field approximation 

The dynamical Hartree-Fock-Bogoliubov (HFB) l29l PiUl 
mean-field equations are recovered by requiring the propaga- 
tion kernels for (f> a and G TabM b to be local in the time and 
space variables. For our model defined by the classical ac- 
tion OTI ) the bare three-point function >5^[0] is given by 
Eq. d57l i and is therefore proportional to <j>. Thus, the HFB 
dynamic equation for the field expectation value cf> is obtained 
by choosing 



and inserting this and T^ 4 ^ hcd [0] = S^J d , see Eq. d40b . into 
Eq. (08]), giving 

iDab( x a,%b)<f>b - -^r- r 4>b(x a )(j>b{Xa)(t)a{Xa) 
C(t a ),b *Al 

= TTTr [Gt a ,bb{Xa,X a )4) a {x a ) + G ta ,ab(x a , X a )4>b{x a ) 



bcl 



- C ( 4 ) 



(69) 



27V 1 



G ta ,ba{Xa,X a )<j>b{Xa)]- 



(70) 



Setting the right-hand side of this equation to zero one recov- 
ers the classical equation of motion. For the free inverse prop- 
agator (f32t this is the Gross-Pitaevskii equation l45ll ; or the 
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non-linear Klein-Gordon equation if D^ 1 is given by Eq. ( l33l ). 
The terms on the right hand side, involving the propagator G, 
are usually referred to as the back action of the excited modes 
onto the condensate within the HFB approximation. 

The HFB evolution equation for G is obtained by neglect- 
ing the term involving 1^ in Eq. d54l ) and inserting the bare 
four-point vertex d40b for 



I 



C{T ab ),C 



iG ^ ac {x ai x c ) - YS°){x a ]<j)) G Tab , cb (x Cl x b ) 

= iSabSc{x a ~ X b ) (71) 



where 



A 

2M 



4>c{Xa) 2 Sab + 24> a (x a )4>b{x a ) 



(72) 



As in the [/(l)-symmetric phase with = the evolution of 
G for x a x b is irrelevant for the evolution on the diagonal, 
x a = x b , such that the equations are entirely local and can be 
rewritten in the well-known form of the time-dependent HFB 
equations, see, e.g., Ref. I26.[29ll. 



3. s-Charmel approximation 

The dynamic equations in s-channel approximation are ob- 
tained, in analogy to the U (l)-symmetric case, by projecting, 
now, the proper 3- and 4-point functions onto the specific ten- 
sor structure, where both functions have at most two different 
external space-time variables x a , xj,. This projection was de- 
fined, for rW, in Eq. ®; for it reads 

r T,alc( X a' X b> X c)[4>] 

= 5t{x a -x b )cj) c (x c )T^ s (x ai x c ) 
+ 5 a c c {x a ~xMxb)T^ s (x ai x b ) 
+ S b c c (x b - x c )Mxa)r^ )s (x a ,Xb). (73) 

The projection implies, in particular, that the four-vertices en- 
tering the full dynamic equation OTb of the propagator, i.e., 
bribed ri 3 L G T , ef if } cd , both have, at most, two dif- 
ferent external space-time variables. Hence, in the flow equa- 
tions for the three- and four-point vertices shown in Fig. [6j 
besides the terms involving 5- and 6-point functions, all terms 
can be discarded with in general more than two external times. 
In the equation for T 1 - 4 ', as before, only the second term on the 
right-hand side is kept, with two 4-point functions. 

As far as the equation for I^ 3 ' is concerned, only the terms 
in the second line in Fig. [6] are kept, with one three, and one 
four-point function, and thus, the fully evolved three-point 
function T T abc is linear in <pd{x) and a sum of terms with 
(d, x) = (a,x a ), (b,Xb), and (c, x c ). As a result, the combi- 

(3) (3) 

nation T T abe G T ^ e f T T j cd is quadratic in (f>, and in the differ- 
ent terms contributing to this vertex, the pair of fields is either 
connected by an internal line, or each is attached to one of the 
external lines, with different space-time variables. 



Choosing the such defined truncation, one derives, as in 
the (J(l)-symmetric case discussed in the previous section, 
coupled integro-differential equations for <f> and G which are 
equivalent to the corresponding set of equations obtained from 
the 2PI effective action in NLO of a 1/7V approximation. 
These equations have been first derived, in Ref. l22ll and can 
be found, e.g., for the case of the model defined in Eq. ( TJTl ), 
in Ref. il. 



VII. CONCLUSIONS 

We have described in detail and extended the approach to 
far-from-equilibrium quantum field dynamics as put forward 
in Ref. J2]. The method is based on a time evolution of the 
system with respect to the flow of a maximum time mark- 
ing the endpoint of the closed time path. In the present work 
we have shown that this approach directly implements the 
underlying causality of the time evolution and has practical 
applications beyond perturbation theory to describe far-from- 
equilibrium quantum dynamics. We have shown in detail that 
in an s-channel truncation the flow equations reduce to dy- 
namic equations as known from a large-A/" expansion of the 
2PI effective action. Although sub-leading terms in the expan- 
sion of the 2PI effective action are suppressed with additional 
powers of 1 /Af they contain bare couplings, and, despite their 
resummation, the approximation seems formally questionable 
for strong couplings if g/Af is large or at large times. An im- 
portant question is to what extent the NLO 1/Af approxima- 
tion is applicable for large interaction strengths. The present 
discussion so far is drawing mainly from benchmark tests for 
special-case systems I27I l4q, l47ll . Here we have shown that 
the self-consistency of such a truncation scheme can be eval- 
uated in a closed form within the present approach by eval- 
uating the flow of the six-point function as induced by the 
four-point function. We emphasize that such a closed self- 
consistency check goes beyond the evaluation of the NNLO 
contribution in a 1/Af expansion. 

Moreover, the equations provided in this work can be eval- 
uated in an iterative manner, beyond the discussed s-channel 
truncation, without the need of resumming specific classes of 
diagrams. This possibility resides on causality which implies 
that the propagation of the integro-differential equations in 
time requires only information computed during the previous 
evolution. We hope to report on this matter in near future. 
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FIG. 6: (color online) Diagrammatic representation of the general 
flow equations for r^ 6c [</)] and r^ 6cd [<£], for a A</> 4 -theory. Open 
circles with a r denote d T R T>a b, solid (blue) lines are r- and, in gen- 
eral, ^-dependent two-point functions G r ,at = i\Tt + R T ]ab ■ AH 
other filled circles denote proper field-dependent n- vertices T^ abcd , 
n — 3,4, 6, using the color scheme of Fig. [4] P implies a sum cor- 
responding to all permutations of its arguments. 
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Appendix A: Flow equations for proper n-point functions 

The flow equations for the one- and two-point proper func- 
tions as derived by functional differentiation of the exact flow 
equation (f30b with respect to the field <f> are given in Eqs. (l35l l 
and ([36}, respectively. In Fig.|6]we provide the flow equations 
for [<j>], for n = 3, 4 and general (f>, in diagrammatic form, 
using the definitions of Fig.@] 



Appendix B: Derivation of the dynamic equations 

In this appendix we derive the dynamic equations from the 
functional flow equations. 



1. Implications from causality 

The specific choice of the sharp cut-off regulator ( TT3T > leads 
to an action functional T T that includes all quantum fluctua- 
tions up to the cut-off time r and reduces to the classical action 
for quantum fluctuations at times greater than r, see Sect. IHTI 

(n) 

For the proper n-point functions T T similar relations hold 
true: 

Property 1 : For a microscopic model with bare vertices that 
are diagonal in time, combined with the sharp cut-off regulator 
( PT5l ). all proper 7i-point functions, with n > 1, obey 

drr^ ai ... a Jx 1 ,...,x n ) = 0, (Bl) 

for r < t niax with i max = max^^o, x n ,o}- This encodes 
the causality of the flow. It follows immediately that 

41...a>l, .••,*«)= itl...^!, • ■ • , In)- (B2) 

Eq. (IB 1 b and (lB2b follow directly from the flow equation < [30b . 
To that end we assume that ( IB2I ) is valid at a given flow time r, 
in particular it holds for r = to- The flow then maintains this 
property. The diagrammatic expansion of the flow equation, 
as e.g. shown in Fig.|6]for n = 3, 4, implies that at least one of 
the external times, t max , is greater than or equal to r for any of 
the diagrams on the right hand side of the flow equation. With 
(IB2I) this particular leg renders the proper function it belongs 
to be equal to its bare counterpart at r = to- Since the bare 
vertices either vanish or are diagonal in the time variables, the 
diagram either vanishes, or at least one time variable of the 
internal propagators attached to it is equal to t max . Since all 
time integrations in the diagrams run from to to, at most r, 
the respective propagators G Tta b(t max , it) vanish. Since the 
other propagators in the diagram, the derivative of the regula- 
tor and the other vertices are bounded or have at most isolated 
singularities (This follows from the fact, that the bare vertices 
have only isolated singularities and that the flow in the cut-off 
time does not give rise to new singularities.) all diagrams on 
the right hand side of the flow equation vanish, which proves 
Eq. ( IB II ). and, by induction and in the limit m —> oo, Eq. ( IB2I ). 

Property 2: For the same kind of microscopic theory and 
regulator as above all proper n-point functions, with n > 1, 
obey Eq. dBlb if r > tj for all j G {1, . . . , n}, and 

^T,l-L...a n (ll ' ■ • ■ ' x n) = ^co,ai...a n (llj • • • ) In)i (B3) 

if t > t max with i max = max{a;i j o, —,x n ,o}. This relation 
states that the proper functions cease to flow once the cut-off 
parameter r exceeds the maximum of the time arguments of 
the correlation function. Eq. dB3b will not be necessary to de- 
rive the dynamic equations from the flow equations such that 
we can use the dynamic equations to derive Eq. ( IB3b . We 
prove this for the example of the truncated equation ( 1331 ) for 
the four-vertex. The time integrals on the right-hand side of 
the dynamic equation 05] ) run from to to the cut-off time t 
and back on the path C. The cut-off times T e f g h appearing 
in the four-vertices on the right-hand side denote the maxi- 
mum of the integration times and are thus smaller or equal to 
t. Expanding Eq. (T55t by iteration, one arrives at a series of 
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diagrams consisting of bare vertices T$ and full propagators 
G Tab ab with cut-off evaluated at the respective maximum of 
time arguments r a b = max(t ai tb). As will be shown in the 
following, in any of these diagrams the contributions from in- 
tegrations of at least one time variable beyond the maximum 
tmax °f external times vanish such that the cut-off parame- 
ter can be set to t = t max which proves Eq. dB3l ). Con- 
sider a diagram with v internal time integrations. The integral 
over the ^-dimensional volume of times in which at least one 
time tj > t max can be divided into integrals over subvolumes 
where the v integration variables are ordered in a particular 
way tpi < tp-2 < ... < tp v where P denotes a permutation 
of the integers 1, v. Any of these subintegrals vanishes sep- 
arately since the respective greatest time t p v is integrated, for 
any value of the second largest time £p(„_i), from ip(„_i) to 
t and, along C, back to £p(„_i). Since the order of integrated 
as well as external times do not change along this integral, 
the time ordering in G does not change and the outward and 
backward contributions along C cancel each other. 



2. G r and G T R T G T for the sharp cut-off 

In this section we prove Eqs. (HTb . ( f42T > which we use in 
Sect. IB 3l to derive the dynamic equations of the proper n-point 
functions. According to Eq. d26l > the propagator G T is the 
inverse of ~i(T^ +R T ). We consider the case that t a , tt, < r. 
As before, R T is the sharp cut-off specified in Eq. (T5\ . Hence 
the integral Jq/^s c R T .ac G TyC b vanishes since R T ^ ac is zero 
for all t c and the propagator is finite. Thus Eq. ( f26b simplifies 
to 



Sc(Xa - X b )5 a b 



r( 2 ) a 



.cb 



(B4) 



C(r),c 



for t a ,tf, < t. Here we also used that both the propagator 
and the proper two-point function are zero for t c > r because 
of Eqs. ( ITTb and dB2L respectively. Therefore the closed time 
path does not contribute for times t c greater than t, and from 
Eq. dB4b one has 



G T ab — i 



(2) 



for T > T ab 



(B5) 



Combining Eq. ( IB5b with Eq. ( fTTb we arrive at Eq. (HTb . 



To derive the identity for (G T [d T R T ]G T ) a b we differentiate 
Eq. ( |26b with respect to r, convolve with the propagator from 
the left. Using again Eq. d26l l we obtain 



/ (2) 
id T G T ,ab = — I G T:ac (d T T T cd )G T , 

JC(oo),cd 

— / G Tyac (d r R T ^ c d)G. 

JC(oo),cd 



ill, 



4b- 



(B6) 



This relation holds for arbitrary regulator functions R T . 

Now, inserting Eq. (|4TT > into Eq. jB6i and substracting 
Eq. (HTT ) differentiated with respect to r, one arrives at 
Eq. (03, where one also uses the consequences of causality, 
Eq. (IB lb . as derived in the previous section. 

3. The dynamic equation for F^ 4 ' 

In the following we present the derivation of the dynamic 
equation for the four-point vertex from the respective flow 
equation given in Eq. d381 l. for the sharp cut-off defined in 
Eq. (TT3T >. To this end we integrate Eq. ( f38l > over the cut-off 
time r from t to t and obtain, after inserting relations ( l4lT> 
and (|42| | for G and GRG, respectively, 



T.abcd 



C(oa),efgh J t a 



dr T 



(4) 

r.abef 
-1 



,(2) 



./'.'/ 



x e{r-t f )9(r-t a )Y^ cdgh 

x [8(t - t h )8{r - t B ) + 9(r - t )6{r - t h )\ 

+ P(a,b,c,d). (B7) 



We use that T 



(4) 

T,abcd 



is symmetric under interchange of the 



a,b,c,d and therefore the integral is invariant under the per- 
mutations (a,b,c,d) — > (b, a, c, d), (a,b,c,d) —> (b,a,d,c), 
and (a, b, c, d) —> (a, b, d, c), of space-time arguments as well 
as field indices. Due to the cyclic property of the trace the 
permutation (a, b, c, d) — > (c,d,a,b) is identical to the re- 
naming of the integration variables (e, /, g, h) — > (g, h, e, /). 
The same is true for the permutations (c, b, a, d) — > (a, d, c, b) 
and (d, b, c, a) — > (c, a, b, d). Using this and performing the 
integration over r we get 
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r (2)" 
. 


-l 
fg 


p(4) 
t f ,cdgh 


r (2)" 


1 6(t 

he 


~ tfMtf 


- te)0(tf 


- tg)0(tf 


-th) 




+ r (4) 
' 1 t g ,abef 


r (2)" 
'» . 


-l 
fg 


r (4) 
t g ,cdgh 


r (2)" 
l tg _ 


0(t 

he 


-t g )0(t g 


- t e )9(tg 


- t f )0(t g - 


-t h ) 




+ r (4) 

T 1 t h .abef 


r (2)" 


-l 
fg 


r (4) 
th,cdgh 


r (2) 


1 6(t 
he 


- t h )9{t h 


- U)9{t h 


- t f )0(t h 


-tg) 



+ (c,b,a,d) + (d,b,c,a), (B8) 
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where the last line denotes two different permutations 
of (a, b, c, d) in the integral. We define T e f gb = 
m&x(t e ,tf,t g ,th) as the maximum of the integration times 
and simplify the above expression using that the sum of prod- 
ucts of ^-functions yields the integration volume J c ^ e ^ gh , 



I =/ 

JC(t),a,...a„ JC( 



Y^e{t-t az ) 6{t a -t aj ). 

IC{t),a\...a n JC(oo),a 1 ...a n i=1 j—ij^i 

(B9) 

Moreover, we use Eqs. d4TT > and (TPTl i to write the integral in 
terms of the propagator G and thus arrive at the form of the 
dynamic equation given in Eq. 



The dynamic equations for r^ 3 ' and ri 4 ^ for i 



The flow equations in the symmetry-broken phase, <f> a 7^ 
0, keeping at most the proper four-point function, are shown 
in Figs. 2] and |6j where the diagrams containing five- or six- 
point functions are omitted. The derivation of the dynamic 
equations for the three- and four-point functions is performed 
as explained in the previous section. We obtain 



,(3) 



I f r (3) r r (3) r r (3) r 

a I 1 r' ,ade T "f' e f r' ,bfg { ~ rT gh,gh L T i chi^nd^d 

D JC(t),defghi 

i 

4 JC(t),efgh 



^rla b efGr fgJg r^\ cgh G Th ^ he + P( ffl> b, c) , 



,(3) 



(BIO) 



r,abcd 



C{t),efghijrs 



p(3) s~i T-i(3) f~< p(3) fi p(3) r-i 

1 t" ,aij^T jB ,]e L T " ,b&f T Ss>*9 t" ,cgh { - r Th r M l r",drs LT T 



r (4) r p(3) r p(3) r 

4 Jc(t), e fghij ' 



I 

Tg 



^labefGr fg JgKl[ cdgh G Th ^ e + P(a, b, C, d) . 



,(4) 



C(t),efgh 



(Bll) 



Here, we use the abbreviations r* = T e fgh, t' — Tdefghi, and 
t" = T e j g hij rs . The truncations made in equations JBlOt and 
( IB 1 U are the omission of the five-point vertex in the dynamic 

(3) 

equation for If and the omission of the five- and six-point 



(4) 



vertices in the dynamic equation for TV 



Appendix C: From the s-channel to the 2PI NLO 1/M 
approximations 

In this appendix we show that Eqs. d65t to (l68l > obtained 
from an s-channel truncation of the dynamic equations de- 
rived in our flow-equation approach are equivalent to the dy- 
namic equations derived from the 2PI effective action in NLO 
of a 1 jM expansion. The proof does not rely on an expansion 
of the integral equations in a perturbative series of Feynman 
diagrams. The steps to prove the equivalence of the s-channel 
and the 2PI NLO 1/M equations are the following: 

Step 1 : We write the integral equation (l65l ) for the s- 
channel vertex in terms of the s-channel coupling function 
(l68l l and split the integral into three parts: The first part con- 
tains no coupling function I T , the second part includes one 
coupling function and one bare vertex, and the third part em- 
bodies two s-channel coupling functions. 



Making use of relations 



and T { t f s (x a ,x b ) 



— (X/J\f)5c(x a — x b ) we rewrite the integral equation 

itM = w L ioo) , cd 

x {0(t - t c )6(t c - t d ) + 6(t - t d )6(t d - t c )} 

x [iI Tcd (x a ,x c ) - S c (x a - x c )]G Ted!cd (x c ,Xd) 

x [U Tcd (x d , Xb) ~ S c (x d - x b )]G Tcd . dc (x dl x c ). (CI) 

We introduced the ^-functions to extend the contour integral 
to infinity in order to reduce the number of different integrals. 
The above integral is now split into three terms: I t = if 5 + 



51 



Til 



removing delta functions by integration, where 



lf\x ai x b )= 0{t-t a )0{t~t b )U Tab {x a ,x b ), (C2) 
if 1 (x a ,Xb) = —i I |/ Tc6 (x a , x c )U Tcb (x c , xb) 



JC(oo),c 

6{t - t c )6(t c - t b )9(t c - t a ) 

f 6{t - t b )6(t b - t c )6(t b - t a 
~h n rac {x a j Xc)I Tca {x c , Xb) 



0(t - t a )6(t a - t c )9(t a - t b ) 
+ 6(t-t c )6(t c -t a )e(t c -t b j\}, (C3) 
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^ T cd ) *^c) n Tcd {x c , X(i) 
C(oo).cd 

x I Tcd (x d ,x b ) 



0(t - t d )9(t d - t c )9(t d - t b )8(t d - t a ) 
+ 9{t - t c )8(t c - t d )8(t c - t b )6{t c - t a ) 
with the basic loop 



(C4) 



A 

W 



G Tab ,cd( x a: x b)G Tab ,cd\ x ai x b)- (C5) 



To obtain the expressions ( IC3b and ( IC4b we renamed integra- 
tion variables to combine all terms into a single integral, and 
we used the property 



I t (x a , x b ) = for all t < r ab 



(C6) 



of the s-channel coupling function to add theta functions as 
follows: 



-^T c d { x a, x c) 



I Tcd (x a ,x c )9(t c - t a ) for t c > t d 
I Tcd {x a ,x c )9(t d -t a ) fort d >t c 



(CI) 

The property ( |C6t follows from the fact that all proper 
?i-point vertices are equal to the respective bare vertex as 
long as one time argument is greater than the cut-off time, see 
Eq. dB2|. 

Step 2: The propagator, i.e., the connected time-ordered 
Green function G Tab . ab (x a ,x b ) = G ab (x a ,x b ) can be de- 
composed into the statistical and spectral correlation func- 
tions F Tab ^ ab (x a , x b ) = F ab {x a ,x b ) and p Tab ,ab{ x a, x b) = 
Pab(x a ,x b ), respectively, 

1 

G ab (x a ,x b ) = F ab (x a ,x b ) - -sgn c (t a - t b )p ab (x ai x b ). 

(C8) 

This decomposition shifts the discontinuity due to the time- 
ordering of the noncommuting field operators into the signum 
function and thus yields two in general independent com- 
ponents, the statistical correlation function defined, for 



bosons, in terms of the anticommutator, F ab (x a , x b ) — 
({&a( x a)> ®b{xb)})c/2 and the spectral function defined 
in terms of the commutator of fields, p a b(x a ,x b ) — 
i([Q a (x a ),$ b (x b )]). The spectral function is related to the 
retarded propagator Gf b (x a ,x b ) = 9(t a - t b )p ab (x a ,x b ) 
as well as to the advanced Green function G^ b (x a ,x b ) = 
-6(t b - t a )p ab (x ai x b ). 

We assume that we can analogously decompose the contour 
s-channel coupling function into two functions in real time by 



I t (x a ,x b ) = lf(x a ,x b ) - -I?(x a ,x b ) sgn c (t a -t b ). (C9) 

It(x a ,x b ) is symmetric under the interchange of x a and x b 
such that lf(x a ,x b ) is symmetric and I p (x a ,x b ) antisym- 
metric under (x a <-> x b ). For the contribution lf s Eq. JC9b is 
obvious from the decomposition of II t , 

Yl^ ab (x a ,x b ) = ^ ^ ab (x a ,x b ) - ^p 2 Tab (x a ,x b ) 

(CIO) 

U r ab ( x a, x b) = j^F TabiCd (x a ,X b )p Tab!cd (Xa,X b ). (Cll) 

For if 1 we show this in Step 3 below, and if 1 vanishes iden- 
tically as proven below. 

We furthermore need that for any x a and x b one has 



Fr( x & 5 x b) 



0. 



(C12) 



For the spectral (p) part of lf s this follows from the relation 

I t ' P (x a ,X b ) = (\/J\f)F TabiCd (x a ,X b )p Tab , cd (x a ,X b ): 

For the relativistic nonlinear Klein-Gordon model 
Pt,cd( x a, x b)\t a =t b = 0. For the Gross-Pitaevskii model 
F is symmetric and p antisymmetric under the exchange 
of field indices, such that the trace over the product F ■ p 
vanishes. For if 1 relation dC12l ) is proven in Step 3. 

In order to prove that // vanishes identically we insert 
Eq. (|C9]) into Eq. d) and use Eqs. ( tCTOl and dCTll to ob- 
tain 



I t ' (x a ,x b ) - -I t ' p (x a ,x b ) sgn c (i a - t b ) = - 



C(oc),cd 



{l? cd (xa,x c ) - l -I p Tc Jx a ,x c ) sgn c (t a - i c )j 



K? cd ( x c x d) - 2 II r C£i ( :z; c^d)sgn c (i c - i d )J yi^ ad {x d ,x b ) - -I p cd (x d ,x b ) sgn c (t d - t b ) 
6{t - t d )9(t d - t c )6(t d - t b )9(t d - t a ) + 9(t - t c )9(t c - t d )9(t c - t b )6(t c - t a ) 

I 



(C13) 



In this expression only the signum functions depend on rewrite the time integrals along the CTP into integrals from 
whether a particular integration time is on the forward or back- to to t. Eliminating the signum functions only terms with 
ward branch of the CTP. The theta-function terms allow to at least two p-components are left. Eq. (IC12I ). the relation 
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W Tah (x a , x b )\ ta=tb = 0, and the fact that 6(t a - t b )9(t b - t a ) 
vanishes except for t a = tb where it is finite then allow to 
show that the integrals over each of the remaining terms in the 
integrand vanish. 



J 



Step 3 : Next we consider the contribution if 1 defined in 
Eq. ( IC3b - Inserting Eq. ( IC9b into Eq. ( IC3b we find its decom- 
positions into contributions lf I F and I t I,p , confirming the 
decomposition 



I t I,F (x a ,Xb) = 0(t - T ah ) 



min(t a ,t(,) 



to 

max(t a ,£(,) 



(l.C 



.\e{t b - t a )I p bc {x ai x c )Ii F cb {x Cl x b ) - 6{t a - t b )IL F ac (x a ,x c )I p ea {x c ,x b )} 



d.r 



.\e{t b - t a )I F c {x a , x c )ILP cb {x c , x b ) - 6{t a - t b )Ii p Tac (x a ,x c )I F a {x c , x b )} 



, (C14) 



lf I,p (x a ,Xb) = 9(t - r ab ) 0{t b -t a ) I dx c I p (x a ,x c )U p (x c ,x b ) - 0(t a - t b ) / dxjl p (x a ,x c )I p ' (x c ,x b ) 



(C15) 



r 



These equations are explicitly symmetric under the exchange 
of x a and Xb- Moreover, we see that Eq. ( IC12I ) is valid for 



as the integral limits are equal in this case. 



T SI,p 

In Eqs. (tC2b . (IC14I ). and ( IC15I ) the theta functions involving 
t imply that the cut-off parameter t of I t (x a , x b ) can be set to 
r ab . Moreover, due to the theta functions in the integrands in 
Eqs. ( IC14l i and ( IC15I ). one can also set the cut-off parameters 
of the / functions appearing in the integrands to the maximum 
of their respective time arguments. Since the cut-off param- 
eters r of any of the functions I T , F T , and p T in the above 
definitions of I F ' P are now fixed to the maximum of their re- 
spective time arguments, they can be neglected. This allows 
to write the integral equations for I^>P{x a , x b ) = I F ' p (x a , x b ) 
as 



and 

Ip(x a ,X b ) = W(x a ,X b ) 



dx c I p (x a ,x c )U p (x c , X b ). 

(C17) 



To arrive at this form one uses that the integrals are 
invariant under the exchange I p (x a , x c )H F (x c , x b ) = 
U F (x a ,x c )I p (x c ,x b ) etc. The integral equations ( IC16I ) and 
( IC17I ) are identical to the equations determining the coupling 
functions I F:P entering the self-energy derived from the 2PI 
effective action in next-to-leading order of a 1 /N expansion, 
see, e.g., Refs. JUIH. 



i F (x a ,x b ) = n F (x Q , x b ) 



dx c I P (x a ,X c ) II F (x c , X b ) 



dx c I F (x a ,x c )U p (x c ,x b ), (C16) 
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